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An integrated system has been developed for discovering potent inhibitors of severe acute respiratory syndrome coronavirus 
3C-like protease (SARS-CoV 3CL pro ) by virtual screening correlating with surface plasmon resonance (SPR) and fluores¬ 
cence resonance energy transfer (FRET) technologies-based assays. The authors screened 81,287 small molecular com¬ 
pounds against SPECS database by virtual screening; 256 compounds were subsequently selected for biological evaluation. 
Through SPR technology-based assay, 52 from these 256 compounds were discovered to show binding to SARS-CoV 3CL pro . 
The enzymatic inhibition activities of these 52 SARS-CoV 3CL pro binders were further applied to FRET-based assay, and 
IC 50 values were determined. Based on this integrated assay platform, 8 new SARS-CoV 3CL pr0 inhibitors were discovered. 
The fact that the obtained IC 50 values for the inhibitors are in good accordance with the discovered dissociation equilibrium 
constants (K D s) assayed by SPR implied the reliability of this platform. Our current work is hoped to supply a powerful 
approach in the discovery of potent SARS-CoV 3CL pro inhibitors, and the determined inhibitors could be used as possible 
lead compounds for further research. (Journal of Biomolecular Screening 2006:915-921) 
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T he last outbreak of the severe acute respiratory 
SYNDROME (SARS) epidemic in 2003 has led to thousands 
of lethally infected patients and hundreds of deaths, 1 and the seri¬ 
ous aftereffect has caused major medical and economic con¬ 
cerns. It has been reported that SARS coronavirus (SARS-CoV) 
is responsible for SARS infection, 2,3 and SARS-CoV 3C-like 
proteinase (SARS-CoV 3CL pro ) is an attractive target for the dis¬ 
covery of anti-SARS agents for its functional importance in 
the viral life cycle. To date, varied kinds of SARS-CoV 3CL pro 
inhibitors were discovered, 4 ' 6 although no effective SARS-CoV 
3CL pro inhibitor has yet been reported to treat SARS. Therefore, 
it is still a great challenge to explore new chemical classes of 


'Drug Discovery and Design Center, State Key Laboratory of Drug Research, 
Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, 
China. 

2 School of Pharmacy, East China University of Science and Technology, 

Shanghai, China. 

Received Mar 28, 2006, and in revised form Jul 7, 2006. Accepted for publica¬ 

tion Jul 18, 2006. 

Journal of Biomolecular Screening 11(8); 2006 

DOI: 10.1177/1087057106293295 


SARS-CoV 3CL pro inhibitors that can be possibly used in 
anti-SARS research. 

As has been used in the current drug-discovery process, by 
using the high-performance computation technique to search the 
large chemical compound databases for the identification of pos¬ 
sible drug candidates, the virtual screening approach is a tech¬ 
nology that is based on the 3D structure of the target protein. It 
involves the rapid fitting of the chemical library members into 
the active sites of 3D protein structures and is critical to distin¬ 
guish active from inactive substances at the primary screening 
stage. To date, there have been successful cases of virtual screen¬ 
ing used in the discovery of lead compound candidates. For 
example, we recently reported that by using virtual screening 
technology with other methods, several new, small molecular 
specific cyclophilin A inhibitors were discovered 7 ; through vir¬ 
tual screening with surface plasmon resonance (SPR) technology 
against a database containing structural information of more than 
8000 available drugs, Cinanserin (SQ 10,643), a well-characterized 
serotonin antagonist that had undergone preliminary clinical test¬ 
ing in humans in the 1960s, was discovered to be the SARS-CoV 
3CL pro inhibitor. 5 

Recently, SPR technology has been recognized as a powerful 
tool in monitoring receptor-ligand interactions with advantages of 
no-labeling, real-time, noninvasive measurements and low sample 
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consumption. 8 It has been widely used in hit determination, 7 ' 9 
and more and more successful cases have been published. 1011 
Fluorescence resonance energy transfer (FRET) technology-based 
assay belongs to another potent system applied in lead compound 
discovery. As a simple and sensitive approach, it has been effec¬ 
tively used in the detection of proteolytic activity and enzyme 
inhibitor evaluation. 4-6 Usually, FRET-based assay involves an 
intramolecularly quenched fluorogenic substrate that contains 
a donor and an acceptor chromophore. In our previous work, 
we constructed an efficient fluorogenic substrate-based platform 
for screening and evaluation of SARS-CoV 3CL pro inhibitors, 12 
with which several potent SARS-CoV 3CL pro inhibitors were 
determined 5 and some related enzymatic features of SARS-CoV 
3CL pro were studied. 13 

In the current work, considering the distinct advantages of 
these 3 above-mentioned methods, and according to the obtained 
successful experiences in enzyme inhibitor discovery, 5 we thus 
developed an integrated strategy for exploring novel SARS-CoV 
3CL pro inhibitors. Our work is expected to supply a sound 
approach in the discovery of new potent SARS-CoV 3CL pro 
inhibitors, and the determined inhibitors could be used as possi¬ 
ble lead compounds for further research. 

MATERIALS AND METHODS 

Materials 

The fluorogenic substrate Dabcyl-KNSTLQSGLRKE-Edans 
was synthesized by Shanghai Sangon Biological Engineering & 
Technology and Service Co. Ltd. (Shanghai, China). The chelat¬ 
ing column and protein molecular weight marker for sodium 
dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) 
use were from Amersham Pharmacia Biosciences (Uppsala, 
Sweden). All other chemicals were of reagent grade or ultrapure 
quality and purchased from Sigma (St Louis, MO). Cinanserin 
was synthesized according to the published method. 5 The tested 
256 compounds from virtual screening were purchased from 
SPECS Company (http://www.specs.net). 

SARS-CoV 3CL pro preparation and activity assay 

SARS-CoV 3CL pro was prepared according to the published 
procedure. 14 The enzymatic activity of SARS-CoV 3CL pro was 
measured by cleavage of the fluorogenic substrate Dabcyl- 
KNSTLQSGLRKE-Edans based on the reported protocols. 5,12 
Briefly, the fluorescence intensity was monitored on a GENios 
microplate reader (Tecan, Mannedorf, Switzerland), and the 
enzyme activity was determined based on the increase in fluo¬ 
rescence upon continuous monitoring of the reactions in 
96-well black microplates (BMG LABTECH, Offenburg, 
Germany) using wavelengths of 340 nm and 488 nm for exci¬ 
tation and emission, respectively. 


Virtual screening 

Virtual screening was performed on a 64-processor SGI 
Origin 3800 supercomputer. The virtual screening protocol 
adopted in this work was similar to that used in the previous pub¬ 
lication. 5 Briefly, the 3D model of SARS-CoV 3CL pi ° was per¬ 
formed as the target for screening the SPECS molecular database 
containing 81,287 molecules by a docking approach 15 and the 
program DOCK 4.0. 16 Residues within a radius of 6 A around the 
catalytic center (His41 and Cysl45) were used to construct 
the scoring grids for the docking screen. During the docking 
calculations, Kollman-all-atom charges 17 and Geisterger-Hiickel 
charges 1819 were assigned to the protein and the small molecules, 
respectively. Conformational flexibility of the compounds from 
the database was considered in the docking search. In DOCK 
simulation, the ligand-receptor binding energy was approxi¬ 
mated by the sum of the van der Waals and electrostatic interac¬ 
tion energies. After initial evaluation of orientation and scoring, 
a grid-based rigid body minimization was carried out for the lig¬ 
and to locate the nearest local energy minimum within the recep¬ 
tor binding site. Position and conformation of each docked 
molecule were optimized using the single anchor search and tor¬ 
sion minimization method of DOCK 4.0. The top 2000 mole¬ 
cules with the highest DOCK score were further subjected to 
CScore 20 and Druglike 21 calculations. Finally, with visual inspec¬ 
tion, 256 compounds were selected and purchased. 

SPR technology-based binding assay and kinetic study 

To inspect the possible active compounds with binding affini¬ 
ties to SARS-3CL pro , SPR technology based Biacore 3000 biosen¬ 
sor was used (Biacore AB, Uppsala, Sweden). SPR assay was 
performed at 25 °C, and SARS-CoV 3CL pro purified at 3 mg/mL 
was diluted to a final concentration of 25 pg/mL in 10 mM 
sodium acetate buffer (pH 4.3) and immobilized to the CM5 
sensor-chip by using the standard primary amine coupling 
method. 22 HBS-EP running buffer (10 mM HEPES, 150 mM 
NaCl, 3.4 mM EDTA, and 0.005% [v/v] surfactant P 20 , pH 7.4) 
was used as the running buffer during immobilization. The sur¬ 
face was activated by injecting a 1:1 mixture containing 0.2 M 
l-Ethyl-3-(3-dimethylaminopropyl) carbodiimide hydrochloride 
(EDC) and 50 mM N-hydroxysuccinimide (NHS) for 7 min at 
a flow rate of 20 pL/min. SARS-CoV 3CL pro was injected and 
the surface was then blocked by injecting 1 M ethanolamine at pH 
8.5 for 7 min. The final immobilization level was about 4000 res¬ 
onance units. 

All the screening assays were performed over the unmodi¬ 
fied dextran surface and SARS-CoV 3CL pro surface. During the 
SARS-CoV 3CL pro binder assay, the sample concentrations 
were set at 1 pM and 10 pM. For sample dilution, the pH and 
the concentration of both DMSO and buffer substances in sam¬ 
ples and running buffer were carefully matched. Each sample 
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assay consisted of a 60-s buffer injection, a 60-s injection, and 
a 120-s undisturbed dissociation phase. The published SARS- 
CoV 3CL pro inhibitor Cinanserin 5 was used as a positive control 
to confirm the binding activity of SARS-CoV 3CL pro , and the 
blank injection was used to check the carryover effects. The 
signal was adjusted for nonspecific binding of the samples to 
dextran matrix by subtracting the signal in the reference chan¬ 
nel from the signal in the active channel. The binding hits could 
be directly determined by the signal difference between the 2 
set concentrations of the samples. 

To further quantitatively investigate the kinetic binding fea¬ 
tures of the determined SARS-CoV 3CL pro inhibitors by FRET- 
based assay, SPR technology-based Biacore 3000 was also 
used. During the assay, samples were diluted to 8 different con¬ 
centrations and injected at a constant flow of 30 pL/min at 
25 °C. Sensorgrams were processed by automatically subtract¬ 
ing for nonspecific bulk refractive index effects. The kinetic 
parameters were then analyzed using a global data analysis 
program (BIAevaluation 3.1 software, Biacore AB). 

SARS-CoV 3CL pro inhibition assay 

SARS-CoV 3CL pro inhibition assay was performed according 
to the previously described method. 5 The inhibition of SARS- 
CoV 3CL pro activity was assayed by preincubating the enzyme 
(1 pM) with the test compounds (50 pM) in 100 mM NaCl, 
20 mM phosphate buffer (pH 7.4) at 4 °C for 2 hours. The reac¬ 
tion was initiated by the addition of 10 pM of fluorescent sub¬ 
strate (Dabcyl-KNSTLQSGLRKE-Edans), and the change of the 
fluorescent intensity was monitored on a GENios microplate 
reader (Tecan) using an excitation wavelength of 340 nm and an 
emission wavelength of 488 nm. Control reactions were carried 
out using the same reaction mixture without the test compound 
or the enzyme. Concentration of inhibitor producing 50% inhibi¬ 
tion (IC 50 ) calculation was carried out in the absence or presence 
of varying concentrations of the inhibitors by using Origin 7.0 
software (OriginLab Corp., Northampton, MA). 

RESULTS AND DISCUSSION 

Virtual screening and docking results 

The DOCK program was employed as the first step in pre¬ 
liminarily screening the potential binding molecules against 
SARS-CoV 3CL pro based on the substrate-binding pocket formed 
by residues within a radius of 6 A around the catalytic center 
(His41 and Cysl45) of the model structure. 23 As a result, 256 
compounds from the SPECS database (containing 81,287 com¬ 
pounds) were thus selected for subsequent biological assays 
according to their DOCK score, Cscore, Druglike score, and 
finally with visual inspection. The results obviously implied that 
the performed virtual screening has supplied potent compound 





FIG. 1. Structures of the 8 compounds screened in this study. 


information, although it might suffer from false positives in view 
of some factors like inaccurate scoring functions, crude confor¬ 
mational statement of target protein and some solvent-related 
terms, 15 or false negatives due to the neglect of allosteric effect as 
carried out in our assay. Shown in Figure 1 are the structures 
of the 8 FRET-determined SARS-CoV 3CL pro inhibitors from 
virtual screening. These compounds show good druglikeness 
according to our in-house druglike filter. 21 Furthermore, 5 of the 
8 compounds pass the rule of 5 based on the results shown in 
Table 1. 

The docking results (Fig. 2) demonstrated that the discovered 
inhibitors can fit well to the active site of SARS-CoV 3CL pro . 
Figure 3 shows the detailed binding information of DC06025 1 to 
SARS-CoV 3CL pro as a typical example. Moreover, by using 
DC060251 as a query structure to search MDL Drug Data Report 
(MDDR), Comprehensive Medicinal Chemistry (CMC), and 
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Table 1. Summary of the Properties of the 
8 SARS-CoV 3CL pro Inhibitors 


Compound 

No. 

Molecular 

Weight 

No. 

Hydrogen 
Bond Donors 

No. 

Hydrogen 
Bond Donors 

CLogP 

Rule of 
5 (Pass 
or Fail) 

DC060015 

558.4 

2 

7 

6.5484 

Fail 

DC060087 

487.5 

1 

9 

2.9486 

Pass 

DC060159 

519.4 

4 

7 

5.6938 

Fail 

DC060170 

541.3 

2 

9 

4.551 

Pass 

DC060180 

522.1 

2 

7 

4.507 

Pass 

DC060245 

536.6 

1 

7 

7.056 

Fail 

DC060251 

470.5 

2 

7 

5.7126 

Pass 

DC060256 

420.3 

1 

5 

3.2525 

Pass 



FIG. 2. The binding model of the 8 inhibitors against SARS-CoV 
3CL pro in the active site. DC060251 is rendered in sticks, whereas others 
are shown in lines. 


China Natural Products Database (CNPD), several known drugs 
or preclinical compounds were also identified (Fig. 4), whose 
structures share 90% similarity to that of DC060251. For example, 
compound A (3,3’,3”,3”’-/?-Xylene-a,a,a’,a’-tetrayltetrakis(4- 
hydroxy-2H-l-benzopyran-2-one)) from MDDR was a known 
HIV-1 protease inhibitor with IC 50 of 1.7 (iM and a mean EC 50 of 
11.5 |iM in an XTT-based anti-HIV assay. 24 Compound B (Ethyl 
Biscoumacetate) from CMC was found in 1943 with anticoagulant 
effect, 25 and compound C (Gerberinol 1) from CNPD was the 
main component extracted from a kind of plant gerbera with many 
biological activities. 26 These results thus implied that compound 
DC060251 might be of great interest for further study. 

SPR technology-based binder identification 

As a 2nd key process, SARS-CoV 3CL pro binder iden¬ 
tification from the virtual screening results was based on SPR 



FIG. 3. Schematic representation of SARS-CoV 3CL pro interaction 
with DC060251. DC060251 and protein residues are rendered in ball- 
and-stick and sticks, respectively. 


technology-based Biacore 3000 biosensor (Biacore AB). 
SPR technology has been widely used to investigate protein- 
macromolecule 27,28 and protein-small molecular compound 5-7 
interactions for its obvious advantages such as label-free, sen¬ 
sitive, real-time, noninvasive measurements, low sample con¬ 
sumption, and high throughput. Recently, this technology has 
also been applied in the discovery of inhibitors against varied 
enzymes, such as HIV-1 protease, 9 human cyclophilin A, 7 and 
human 5-lipoxygenase. 10 In our previous reports, we ever suc¬ 
cessfully used SPR technology for discovering SARS-CoV 
3CL pro small molecular inhibitors 5 and performing SARS-CoV 
3CL pro -related research. 13 

In the current work, the SPR-based results showed that of the 
selected 256 compounds from the virtual screening, 52 com¬ 
pounds were identified as SARS-CoV 3CL pro binders. SPR- 
based technique has supplied a rapid approach in the discovery 
of SARS-CoV 3CL pro binder. However, it is noticed that in some 
cases some hits obtained by SPR assay are probably nonspecific 
binding molecules, which show no inhibitory activities. To solve 
this problem, the screened hits should be further tested by enzy¬ 
matic inhibitory assays. 

FRET technology-based IC 50 determination 

As demonstrated in the SPR-based SARS-CoV 3CL pr0 binder 
assay, 52 SARS-CoV 3CL pro binders were discovered. To further 
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FIG. 4. Structures of 3 compounds (A. B, C) from MDL Drug Data 
Report (MDDR), Comprehensive Medicinal Chemistry (CMC), and 
China Natural Products Database (CNPD) databases. DC060251 
shows 90% structural similarity to them. 


investigate their possible inhibitory activities against SARS-CoV 
3CL pro , the inhibitory effects of the selected 52 binders were 
tested by the FRET method. For this assay, the internally 
quenched fluorogenic substrate Dabcyl-KNSTLQSGLRKE- 
Edans was applied, whose Michaelis constant (Km) against 
SARS-CoV 3CL pro is 49.38 ± 5.80 jxM (data not shown). This 
value is comparable to those of other recently published fluo¬ 
rogenic peptides. 6,29 

Table 2 lists the yielded IC 50 values of the 8 screened inhibitors 
of SARS-CoV 3CL pro , and Figure 5 shows the dose-response 
curves for 3 typical inhibitors DC060170 (A), DC060180 (B), and 
DC060251 (C). The IC 50 values of these inhibitors are in the range 
of 6.86 to 80.46 |xM. 


SPR-based quantitative inhibitor binding affinity 
determination 

It has been known that although the IC 50 value represents the 
inhibitory activity of an inhibitor to an enzyme, this value is 
subject to the experimental conditions, for example, the struc¬ 
ture of the substrate and the total concentrations of both enzyme 
and substrate. If any of the experimental conditions change, the 
IC 50 value of a given compound will be altered, making IC 50 
from different literatures hard to be compared. 30 In comparison 
with IC 50 , the equilibrium dissociation constant (K D ) for the 
enzyme-inhibitor complex is independent of the experimental 
conditions. SPR-based assay can directly measure K D values 
between the inhibitors and SARS-CoV 3CL pro . To further inves¬ 
tigate the K d values of the inhibitors to SARS-CoV 3CL pro , an 
SPR technology-based Biacore 3000 instrument was also used. 
K d values of the 8 discovered inhibitors were calculated using 


Table 2. Summary of Concentration of Inhibitor 
Producing 50% Inhibition (IC 50 ) Values and 
Equilibrium Dissociation Constants (K D s) 
for SARS-CoV 3CL pro Inhibitors 


Compound 

Molecular 
Weight (Da) 

IC S0 (gM) 

K D (gM) 

DC060015 

558.4 

49.96 ± 4.62 

18.60 ± 1.18 

DC060087 

487.5 

80.46 ± 2.07 

42.50 ±2.10 

DC060159 

519.4 

42.19 ±3.02 

19.14 ±2.03 

DC060170 

541.3 

6.86 ± 0.91 

4.23 ±0.81 

DC060180 

522.1 

9.19 ±2.34 

6.50 ± 0.23 

DC060245 

536.6 

21.49 ±2.97 

23.91 ±5.63 

DC060251 

470.5 

9.56 ± 1.75 

4.78 ± 0.29 

DC060256 

420.3 

56.24 ± 3.48 

49.91 ±3.98 


the 1:1 Langmuir binding model or the steady-state affinity 
model (BIAevaluation 3.1 software, Biacore AB). Figure 5 
(inset) shows the binding dose-response curves for the 3 com¬ 
pounds, DC060170 (A), DC060180 (B), and DC060251 (C). As 
indicated in Table 2, the K D values of the compounds range 
from 4.23 to 49.91 jiM. As demonstrated in Figure 6, the K D 
values are well correlated with IC 50 values by the correlation 
coefficient ( R ) of 0.92. 


CONCLUSIONS 

In this work, we developed an integrated approach for dis¬ 
covering novel SARS-CoV 3CL pro inhibitors. Virtual screening 
method based on the 3D structure of SARS-CoV 3CL pro can 
reliably narrow down the number of potential candidates before 
experimental testing. Although virtual screening might suffer 
from false positives, SPR technology was thus used to deter¬ 
mine the binders of SARS-CoV 3CL pro and further reduce the 
screening scale and time. Because the hits obtained by SPR- 
based assay may be nonspecific binding compounds, FRET 
method is applied for SARS-CoV 3CL pro inhibitor identifica¬ 
tion. Our current strategy for discovering SARS-CoV 3CL pro 
inhibitor has integrated virtual screening and SPR- and FRET- 
based assays, and with this integrated platform, we have suc¬ 
cessfully identified 8 novel potent SARS-CoV 3CL pro inhibitors. 
Although the hit enrichment might not always be higher than a 
random selection, the screened active compounds from virtual 
screening could involve more structural information for virtual 
screening based on the compound database that comprises a 
large number of chemical structures, whereas for random 
screening, the obtained active compounds might contain lower 
structural information because of the possible limitation of the 
compound quantity. The obtained results thereby suggested that 
the inhibitor screening is more effective compared with the tradi¬ 
tional screening when integrated strategy is applied. It is 
believed that such an integrated approach may help speed up the 
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Log(C) M 

FIG. 5. Dose-response inhibition and kinetic assay curves for the 3 
compounds from the screening procedure. IC 50 values were extracted 
from the curves using nonlinear regression analysis, and dissociation 
equilibrium constant (K D ) values were analyzed using a global data 
analysis program of Biacore 3000. A, B, C represent the dose-response 
inhibition curves (Inset: the responding binding curve) of the compounds 
(DC060170, DC060180, DC060251), respectively. 
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IC 50 (nM) 

FIG. 6. Correlation between the concentration of inhibitor produc¬ 
ing 50% inhibition (IC 50 ) and equilibrium dissociation constants (K D ) 
for SARS-CoV 3CL pro inhibitors. IC 50 values were obtained from 
the fluorescence resonance energy transfer (FRET)-based assay, and 
K d values from biosensor measurements. The symbol (•) represents 
the compounds with correlation. 


drug-discovery process. In addition, it is expected that the deter¬ 
mined inhibitors could be used as possible lead compounds for 
further research. 
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